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Abstract 


The X-ray luminosity function of clusters of galaxies is determined at different cosmic 
epochs using data from the Einstein Observatory Extended Medium Sensitivity Survey. The 
sample consists of 67 X-ray selected clusters that have been grouped in three redshift shells. 
Evolution is detected in the X-ray properties of clusters. The present volume density of high 
luminosity clusters is found to be greater than it was in the past. Given the still limited data 
set this result should be regarded as preliminary. It can be interpreted as the consequence 
of either luminosity evolution or modest density evolution. 

Keywords: Clusters of Galaxies: Luminosity Function - Evolution - X-rays 
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I. Introduction 


The study of distant clusters of galaxies provides important information on their forma- 
tion and evolution. Investigations of the X-ray evolution have almost always proceeded by 
making observations of clusters of galaxies selected in the optical (Henry et al. 1982; Henry 
and Lavery, 1984). There is currently little evidence for evolution in X-ray luminosity or 
temperature for distant optically selected clusters. Surprisingly, for the only two distant sys- 
tems studied in detail, the cluster 0016+16 at z = 0.541 (White, Silk and Henry, 1981) and 
the cluster around 3C295 at z = 0.461 (Henry and Henriksen, 1986), the X-ray properties 
were found to be similar to those of nearer rich clusters. However, these apparent similari- 
ties between distant and nearby systems might be primarily due to a selection effect. Since 
distant clusters selected optically are chosen because they are especially rich, these clusters 
may be among the few which have already undergone a considerable amount of dynamical 
evolution. It is almost impossible to avoid or quantify biases in optically selected samples 
because they are chosen by eye. Even investigators who presently make catalogs by scanning 
plates and who select the galaxies and clusters by using rigorous algorithms are faced with 
the problem of contamination by foreground galaxies and stars. 

X-ray selection does not have these biases, even though different selection effects are 
present. There may be a preference for the detection of high surface brightness systems as 
well as clusters with deep potential wells. Since the vast majority of the clusters known today 
have been selected in the optical, it is vital to investigate the properties of a cluster sample 
extracted from, an X-ray survey for a different approach to the understanding of cluster 
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formation and evolution. A Hubble constant H 0 = 50 km s" 1 Mpc \ and a Friedmann 
universe with a deceleration parameter q 0 = 0 is assumed throughout this Letter. 

II. The Sample 

The sample of clusters of galaxies used in this study is extracted from the Einstein 
Extended Medium Sensitivity Survey (EMSS). A detailed description of the survey sources, 
the selection criteria, the detection algorithm, the computation of the X-ray flux and other 
parameters is given in Gioia et al. 1990. We recall here that the EMSS is a flux limited 
sample consisting of 835 sources serendipitously detected in Imaging Proportional Counter 
(IPC) fields at high Galactic latitude, with limiting sensitivities ranging from 5 x 10" 14 
to 3 x 10~ 12 erg cm -2 s -1 in the 0. 3-3.5 keV energy band. A detailed discussion of the 
identifications, as well as presentation of finding charts, spectral and photometric data is in 
preparation and will be presented elsewhere (Stocke et al. 1990, in preparation; Maccacaro et 
al. 1990, in preparation). For the purpose of this study we have considered only sources with 
declination greater than -40° (accessible from Mauna Kea) and flux greater than 1.5 x 10 -13 
erg cm -2 s -1 in a 2/4 x 2/4 detection cell (to reduce the number of still unidentified sources). 
Adopting these criteria the survey contains 733 sources and is 97% identified. There are 93 
sources identified with clusters of galaxies. Since most nearby clusters were observed as a 
target of IPC observations they were not available to be detected serendipitously by the 
EMSS, so this sample is not complete at the low end of the redshift distribution. For this 
study we have chosen to use only those clusters in our sample with a redshift greater than 
0.14. This value roughly corresponds to Abell Distance classes 5-6. Since the majority of 
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the Abell clusters chosen as the target of the IPC observations belong to Distance Class 3 or 
less, we feel comfortable in using the value of 0.14 as a lower limit in redshift. The resulting 
sample contains 67 objects. It is the most numerous sample of distant clusters of galaxies 
extracted from a flux limited survey of “faint” X-ray sources: i.e. the sample is defined 
exclusively by its X-ray properties. The precise knowledge of the area of sky searched for 
X-ray sources and of the limiting sensitivity pertaining to each area allows us to derive the 
cluster X-ray luminosity function. In the redshift range 0.14 < z < 0.20 there are 20 clusters 
of which only 4 are in the Abell catalog. The X-ray luminosities of objects in this shell are 
all greater than 10 44 erg s _l . Since clusters with this luminosity are almost exclusively Abell 
clusters at lower redshift, it is somewhat surprising that we find so few of them. Although 
our clusters are mostly not Abell clusters, it is premature to discuss the implications of these 
results on the completness of the Abell catalog until more optical work is completed on our 
clusters. 

Even though the EMSS is statistically well defined, there are still a number of effects in 
the data which must be taken into account. These effects are absorption by the Milky Way, 
the different redshifts of the sources, the different sky coverage for different flux limits, and 
the correction for lost flux due to the finite source size (the EMSS uses a detection cell of 2/4 
x 2/4). We discuss the corrections we have applied to take into account each of these effects 
on our data. We note here that the EMSS uses the so called M-DETECT algorithm to find 
sources. In this method the background is computed from a global map of the detector so 
that sources are not lost because their extended flux distribution mistakenly increases the 
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apparent background around them (see Gioia et al. 1990 for a detailed discussion). 

The flux from each source has been corrected for absorption using the neutral Hydrogen 
values from the survey of Stark et al. (1984). Most of the sky was observed through a small 
range of N// which results in a negligible bias (see Zamorani et al. 1988, and Maccacaro et 
al. 1988). K-corrections are small for our sample. Assuming a Raymond-Smith spectrum at 
a redshift of 0.5, the correction is always within 15% of unity for temperatures between 2 
and 10 keV (Burg, 1989 private communication) and is less at smaller redshifts. Therefore, 
for simplicity we used K-corrections calculated assuming a power law spectrum with energy 
index of 0.5, which roughly approximates a 6 keV thermal spectrum in our 0.3— 3.5 keV 
energy band. The sky coverage North of —40° declination has been calculated adopting the 
same spectrum and procedures used to determine the flux of the clusters, i.e. using only the 
counts in the detection cell without applying any correction for the point response function 
or the mirror scattering. 

The largest correction is for the effect of the finite size of the X-ray emission in clusters. 
This correction will be used in two places: first to calculate the true luminosity of a cluster 
when only its detection cell flux and redshift are known, and second to determine the amount 
of flux that would appear in the detection cell at different redshifts during the calculation 
of the maximum observable redshift (see Section III) in the derivation of the luminosity 
function. We adopt the (3 model for the cluster surface brightness, that is 

= [1 + (0/0o) 2 ] (3/, - 1/2) 

From Jones and Forman (1984) we adopt (3 = 2/3. Then integrating over the square detection 
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cell between 0 and 9 D we obtain the observed flux: 


f$D f&D 

F ob , = 4 / d0j dO y I{0) 
Jo Jo 


where 9q is the angular half-size of the detection cell (1/2) we obtain: 


F oh , = 27t/ o 0 o 2 /(^) 
00 


where 


n e 0 ] ~ 2 + k 1 (Sl/ol + iy 1 


( 1 ) 


is the fraction of the total flux in the detection cell. We determine 9 0 using the 18 Piccinotti 
et al. (1982) clusters (HEA01-A2 experiment) which have IPC imaging data and which 
are not so large that they extend beyond the IPC ribs. This sample is X-ray selected and 
seems to be the most comparable with our sample. For these 18 clusters, we calculate the 
average fraction of the total flux (as determined by Piccinotti et al.) which is detected by the 
EMSS detection cell if each cluster were at a fiducial redshift. The arbitrary fiducial redshift 
was chosen to be 0.35 which gives D a x9 d = 0.5 Mpc, where D A is the angular diameter 
distance. At this redshift, the average fraction of flux of the Piccinotti et al. objects that 
would be in the detection cell is equal to 0.43 with a large dispersion (a = 0.2). No obvious 
dependence is found between this fraction and the X-ray luminosity of the clusters. From 
(1), D A x 9 0 = 0.37 Mpc. Equation (1) with this value of 9 0 is used to correct the observed 
luminosity to the total luminosity. The redshift dependence of Iq and 6q are: 

7q(1 + ZobaY 


h (z) = 
9 0 (z) = 


( 1+*) 4 

9qD A (z oba ') 

~dazT 
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where z ob , is the redshift of the cluster. We assume that 9q and do not evolve with redshift. 
Performing a similar integral over the detection cell for the cluster at an arbitrary redshift 
gives 


F(z) = F obs 


Dl(z oba ) /(DaWd/DM 


( 2 ) 


Dl(z) f(D A (z obs )9 D /D A 9 0 ) 
where Di is the luminosity distance. This expression is used in the calculation of the 
maximum redshift at which a given object could have been detected. It reduces to the point 
source result in the limit that the size of the detection cell is much larger than that of the 


cluster. 


III. The X-ray Luminosity Function 

In this section we derive and discuss the X-ray Luminosity Function (XLF) for clusters 
of galaxies computed in three different redshift shells defined by z\ ow and z high- A non- 
parametric representation of the XLF is first obtained using the 1/V a method of Avni and 
Bahcall (1980) a generalization of the 1/V maz method (Schmidt 1968) when several complete 
samples are analyzed. For each cluster falling in one of the three ziow—Zhigh shells defined 
later we computed the maximum redshift (z max ) at which the source could have been seen 
taking into account the solid angle observed at different limiting sensitivities. The maximum 
volume in which the cluster could have been detected depends on the redshift of the shell 
under consideration, the luminosity of the cluster and the sky coverage of the EMSS. The 
search volume for a given cluster, V a , is the sum of all volumes lying between the minimum 
redshift z/ ou> of the shell under consideration and the lesser of the maximum shell redshift z high 
or the maximum redshift z mor at which the source could have been seen for each different 
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sensitivity limit. The individual contributions have been binned by log luminosity to create 
the differential XLF N(L) integrated in independent bins 0.3 wide in AlogL. For each bin 
we have: 


N(L) = 


n 


E 


i 

v a ,AL 


where n is the number of objects in that bin. The results are shown in Fig. 1, where 
the three panels give the XLF’s in redshift shells as indicated. There are 20 clusters with 
redshift 0.14 < 2 < 0.20, 26 clusters with redshift 0.20 < z < 0.30 and 21 clusters with 
redshift 0.30 < z < 0.60. The lcr error bars associated with each bin are determined 
from the number of objects contributing to that bin and have been computed using Poisson 
statistics (Regener, 1951). 


We then consider a parametric representation of the luminosity function of the form 

-A- = I<L; t °Mpc- 3 L;t 

dLi \ 4 

where L 44 is the X-ray luminosity in units of 10 44 erg s -1 and K is the normalization coefficient 
expressed in units of Mpc -3 L 44 . The maximum likelihood method (see Murdoch, Crawford 
and Jauncey, 1973, and references therein) has been applied to the unbinned data to estimate 
the best fit slopes which are given in Table 1 with there associated 1 a errors. The fits have 
been computed in each redshift range between the minimum observable luminosity L mtn and 
infinite luminosity. 
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Table 1 


Parametric Representation of the Cluster X-ray Luminosity Function 


z 

a 

K{Mpc~ 3 L ^) 

log L m { n 

0.14-0.20 

2.09 ± 0.20 

(7.19 ±0.58) x 10- 7 

42.90 

0.20 - 0.30 

2.63 ± 0.22 

(10.8 ±1.56) x 10- 7 

43.30 

0.30-0.60 

3.09 ± 0.27 

(12.2 ±4.46) x IQ" 7 

43.80 


The normalization coefficient K has been computed by requiring that the number of 
expected objects equals the number of observed objects. Errors on K have been determined 
by letting a assume the la extremes in each case. A steepening of the slope is observed at 
higher redshifts. This change is best seen in Fig. 2 where the differential XLF’s for clusters 
in the lowest shell (0.14 < z < 0.20) and in the highest shell (0.30 < z < 0.60) are plotted. 
The difference between the two slopes is significant at the 3 a confidence level. 

We note that no cluster with log L r > 45.2 has been detected at low redshift (0.14 < z < 
0.20). Clearly these clusters, which have been detected at higher redshift, could have been 
detected at lower redshift. However the available volume in the low redshift shell is much 
smaller than the volume in the higher redshift shells. Only one cluster is expected in the bin 
centered at log L x = 45.35. Thus the absence of log L r > 45.2 clusters in the low-z shell is 
not significant and is not necessarily indicative of a break in the XLF of low redshift clusters. 
It is of interest to compare our lowest redshift XLF with that of Piccinotti et al. (1982) even 
though their sample was at a lower redshift and in a different energy band. The two slopes 
agree (2.09 vs. 2.15) within the errors. With the same 6 keV thermal spectrum we used 
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previously to compute our fluxes, we find that the Piccinotti norma.liza.tion, converted from 
their 2-10 keV band to our 0.3-3.5 keV, is 2.1 ± 1.6 times smaller than ours. That is we 
essentially agree within the errors. 

IV. Discussion 

Even if our result is significant only at the 3cr level we believe we have the first con- 
vincing evidence for evolution in the X-ray properties of clusters. Note that the slopes are 
independent of the correction applied for the flux lost outside of the detection cell because 
all the sources in each shell are at approximately the same redshift. The XLF which char- 
acterizes high redshift clusters (0.30 < z < 0.60) is significantly steeper than the XLF of 
low redshift clusters (0.14 < z < 0.20). This trend is supported by the intermediate XLF 
(0.20 < z < 0.30). We have assumed that 9 0 does not evolve. If, however, the core radius 
decreases with redshift, as would be expected in an hierarchical scenario, then we would 
make a smaller correction for the flux lost outside the detection cell, would have fewer high 
luminosity clusters, and would find even stronger evolution than we do. 

The luminosity range covered by our data is rather limited and prevents a detailed anal- 
ysis of the shape and kind of cosmological evolution. Furthermore, there are still 19 sources 
unidentified which could modify the results presented in this Letter. However evolution seems 
present only at high luminosities (log L r > 44.7). At lower luminosities (44.2 < logL x < 
44.7) no significant difference exists as a function of redshift. This behavior is suggestive of a 
luminosity dependent evolution such that the volume density of high luminosity clusters (log 
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L x > 44.7) is higher at the present epoch than at epochs corresponding to redshifts of about 
0.5. The volume density of low luminosity clusters has remained unchanged. Presumably 
our power law characterization of the data could be naive and a luminosity function with a 
break analogous to a Schechter function could be more appropriate. In this case, the differ- 
ent slopes observed could be the result of either luminosity evolution , with the break shifted 
to a higher luminosity for the lower redshift clusters, or a modest density evolution. In the 
latter case the apparent change in slope would result because the high redshift sample, with 
intrinsically more luminous clusters in the mean, has more objects drawn from above the 
break relative to the low redshift samples. 

The basic conclusion is that there is good evidence for a difference between the X-ray 
luminosity function at high and low redshift. This difference goes in the opposite sense to 
that anticipated by Kaiser (1986) who predicted density evolution in the sense that there 
would have been a higher volume density of X-ray clusters in the past. By contrast, the 
evidence for evolution we have found, that is fewer clusters in the past, is in the same sense 
as anticipated by Perrenod (1980) (see also Cavaliere and Colafrancesco, 1988). 

After this paper had been submitted we received a preprint from Edge et al. (1990) who 
found a similar result from an independent data set. 

We would like to thank R. Giacconi and R. Burg for many interesting discussions when this 
work was at an early stage of preparation, and G. Zamorani for many helpful suggestions 
and comments. This work has received partial financial support from NASA contract N AS8- 
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Figure Captions 


Fig. 1 - The differential luminosity functions of clusters in the three redshift shells: a) 
0.14 < 2 < 0.20; b) 0.20 < z < 0.30 and c) 0.30 < z < 0.60. 

Fig. 2 - A comparison of the differential luminosity functions for the lowest (open squares) 
and the highest (filled dots) redshift shells. The straight lines represent maximum 
likelihood fits to the individual unbinned data. The dot-dashed line is the fit for the 
lowest redshift shell and the solid line is the fit for the highest redshift shell. 
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